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ABSTRACT 


The  absolute  vorticity  and  angular  momentum  budgets  for 
the  Presidents'  Day  Cyclone,  18-19  February  1979  are  evalu- 
ated employing  the  guasi-Lagrangian  diagnostic  technique  in 
pressure  coordinates.  The  First  GARP  Global  Experiment 
(FGGE)  Level  Ill-b  data  set  prepared  by  the  European  Center 
for  Medium  Range  Weather  Forecasts  is  used  to  define  the 
storm  and   to   calculate   the    budgets. 

Important  differences  are  found  between  the  F3GE  and 
mass  budget  vertical  motion  estimates.  The  budget  results 
illustrate  the  importance  of  inward  eddy  node  lateral  trans- 
port and  advection  at  upper  levels  and  lean  mode  transport 
due  to  convergence  at  lower  levels.  The  vertical  redistri- 
bution plays  a  lesser  role  but  is  more  significant  in  the 
angular  momentum  budget.  The  vertical  vorticity  advection 
and  tilting  terms  are  comparable  in  magnitude.  The  largest 
budget  residuals  occur  during  the  period  of  most  intense 
convective  activity  which  indicates  a  need  for  a  better 
understanding  of  the   vertical   redistribution   processes. 
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I.        INTRODUCTION 

Operational  forecast  models  have  generally  performed 
poorly  in  predicting  the  frequently  observed,  cold-season 
explosive   cyclogenesis      of    maritime  cyclones.  Sanders   and 

Gyakum  (1980)  demonstrated  that  the  National  Meteorological 
Center's  (NMC)  Limited  area  Fine-Mesh  Model  (LFM)  predicts 
only  one-fourth  to  one-third  of  the  deepening  rates  of  rap- 
idly  deepening  maritime   cyclones. 

The  accurate  forecasting  of  these  meteorological  events 
over  the  open  ocean,  as  well  as  in  coastal  regions,  is  of 
critical  importance  to  both  civilian  and  military  interests. 
The  extreme  conditions  produced  by  these  storms  often  result 
in  less  of  life,  destruction  of  property  and  disruption  of 
operations.  The  recent  sinking  of  the  oil  rig  Ocean  Ranger 
(LeMoyne,  1982)  and  the  attendant  loss  of  life  illustrates 
the  magnitude  of  the  destructive  forces  oreated  by  intense 
maritime      cyclones.  Throughout  recorded      history,         these 

storms  have  claimed  numerous  ocean-going  vessels  in  prefer- 
red areas  of  explosive  cyclogenesis  such  as  off  the  east 
coast  of  the  U.S.  These  cyclones  often  strike  with  little 
or  no  warning.  More  accurate  and  timely  forecasting  of 
explosively  developing  maritime  cyclones  would  undoubtedly 
reduce  their   impact  on   human   activity. 

11 


The  failure  of  numerical  forecast  models  in  predicting 
the  intensities  of  vigorous  maritime  cyclones  is  due  to  the 
small  scales  over  which  the  explosive  development  occurs  as 
well  as  the  lack  of  meteorological  observations  over  oceanic 
areas.  This  problem  has  resulted  in  an  inaccurate  represen- 
tation of  the  initial  conditions  in  the  forecast  models. 
The  most  promising  solution  to  this  problem  is  the  incorpo- 
ration of  wind,  temperature,  moisture  and  pressure  data 
obtained  from  satellite  measureients  in  the  operational  pre- 
diction models.  It  appears  however,  that  a  better  under- 
standing of  the  role  of  dynamic  forcing  in  maritime 
cyclogenesis  is  required  if  this  new  data  source  is  to  be 
fully  utilized. 

This  study  constitutes  a  portion  of  a  broad-based 
research  effort  zo  better  understand  the  processes  leading 
to  the  development,  maturation  and  decay  of  extratropical 
maritime  cyclones.  This  improved  understanding  will  point 
toward  ways  to  improve  the  numerical  prediction  of  maritime 
cyclogenesis. 

More  specifically,  this  thesis  consists  of  circulation 
(vorticity)  and  angular  momentum  budget  analyses  of  the 
Presidents1  Day  Cyclone  (18-19  February  1979);  a  subsynoptic 
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scale  event  which  provided  a  well-defined  example  of 
explosive  maritime  cyclogenesis .  This  storm  has  been  the 
subject  of  several  in-depth  analyses  (Bosart,  1981;  Roman, 
1981;  Uccellini,  Kocin  and  Mash,  1981).  It  produced  record 
breaking  snowfalls  (nearly  60  cm)  in  portions  of  Virginia, 
Maryland      and      Delaware.  The        central      surface      pressure 

decreased  by  more  than  27  mb  luring  a  24-hour  period.  The 
NMC  LFM-II  model  underf or ecast  the  decrease  in  the  surface 
pressure  minimum  by  an  average  of  8  to  1 6  mb  during  the 
explosive  cyclogenesis  period  (Bosart,  1981).  Thus,  the 
Presidents1  Day  Cyclone  can  be  included  ia  a  class  of  explo- 
sive deepeners  not  adequately  handled  by  current  numerical 
forecast  models  (Sanders  and  Syakum,  1930) .  This  cyclone 
developed  during  one  of  the  First  GARP  Global  Experiment 
(FGGE)  special  observing  periods.  Consequently,  there  are 
additional  observational  data  and  research-based  objective 
analyses    to  aid    in  the  study. 

Roman  (1981)  completed  a  mass  budget  study  of  the  Presi- 
dents1 Day  Cyclone  based  on  the  FGGE  data  set.  He  computed 
vertical  velocity  values  in  layers  consistent  with  the  lat- 
eral mass  transport  and  time  tendency  terms.  The  vertical 
velocity    fields    in   the   FGGE  data    set    are  used   to   compute  the 
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vertical  transport  terms  in  the  budget  studies  in  this 
thesis. 

Quasi- Lagrangian  diagnostic  (QLD)  procedures  (Wash, 
1978)  are  employed  in  this  study.  The  basic  feature  of  the 
QLD  budget  approach  is  the  computation  of  volume  integrals 
of  atmospheric  properties.  This  volume  is  centered  on  and 
moved  with  the  surface  pressure  minimum  of  the  Presidents1 
Day  Cyclone.  This  approach  incorporates  the  translational 
effects  associated  with  the  cyclone,  as  wall  as  the  effects 
due  to  development  or  decay.  The  lateral  exchanges,  sources 
and  sinks,  and  vertical  redistributions  of  meteorological 
variables  are  primarily  the  result  of  developmental  pro- 
cesses. QLD  budgets  yield  quantitative  evaluations  of  the 
exchanges  of  mass,  circulation,  angular  momentum  and  kinetic 
energy  between  the  budget  volume  and  its  environment. 

The  circulation  about  the  surface  cyclone  increases  dra- 
matically during  cyclogenesis,  and  provides  a  convenient 
measure  of  cyclone  intensity.  The  circulation  budget  is 
used  to  evaluate  the  magnitude  of  the  divergence  and  hori- 
zontal and  vertical  advection  terms  in  the  vorticity  equa- 
tion. It  is  well  known  that  upper-level  vorticity  advection 
contributes  to  cyclogenetic  processes  via  the  divergence  and 


14 


vertical  velocity  fields  (Petterssen,  1955).  This  qualita- 
tive relationship  does  not  howsver,  pro/ide  answers  to  key 
questions  such  as: 

•  What   is  the   magnitude  and   significance  of  low-level 
horizontal  vorticity  advection? 

•  How  does  the  vertical  advection  of  vorticity  contribute 
to  cyclogenesis? 

•  What  is   the  magnitude  of   the   local   generation   of 
vorticity? 

The  quantitative  analysis  of  circulation  using  budget  diag- 
nostics provides  answers  to  these  questions  through  the 
identification  of  lateral  exchanges,  vertical  redistribu- 
tions and  sources  and  sinks. 

An  alternative,  but  closely  related,  measure  of  vortex 
intensity  is  that  of  angular  momentum.  The  concentration  of 
momentum  about  the  storm  axis  is  primarily  the  result  of 
radial  transport.  Johnson  and  Downey  (1975)  show  that  there 
are  no  sources  of  angular  momentum  within  the  moving  budget 
volume.  Large  scale  forcing  of  storm  angular  momentum 
buildup  occurs  through  inward  transport  of  angular  momentum 
at  lower  levels.  One  objective  of  the  QLD  analysis  of  angu- 
lar momentum  is  to  reveal  the  significance  of  the  vertical 
transport  of  this  property  daring  the  development  of  the 
Presidents'  Day  Cyclone. 
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The  budget  analyses  are  dona  in  isobaric  vertical  coor- 
dinates. The  choice  of  these  coordinates  is  made  for  con- 
venience in  comparing  the  budget  analyses  with  standard 
diagnostic  and  prognostic  charts.  These  analyses  can  also 
be  converted  to  sigma  vertical  coordinates  over  level  ter- 
rain to  facilitate  comparisons  with  resuLts  from  numerical 
model  predictions. 

In  summary,  the  purpose  of  this  thesis  is  to  examine  a 
well-defined  maritime  cyclogenasis  event  in  terms  of  circu- 
lation and  angular  momentum  from  the  budget  perspective. 
The  specific  objectives  of  this  study  are: 

•  Quantitatively  evaluate  the  horizontal  transport,  ver- 
tical redistribution  and  sd urces/sinks  of  vorticity  and 
angular  momentum  within  a  budget  volume  centered  on  the 
Presidents1  Day  Cyclone; 


•  Compare  the  vertical  velocity  structure  of  the  Presi- 
dents1 Day  Cyclone  derived  from  the  mass  budget  verti- 
cal transport  (Roman,  1981)  with  that  of  the  F3GE  data 
set. 

A  brief  description  of  the  FGGE  data  set  and  a  synoptic 
overview  of  the  Presidents*  Day  Cyclone  are  presented  in  the 
following  chapter.  Chapter  III  contains  an  abbreviated 
description  of  the  QLD  technique.  A  mora  detailed  explana- 
tion of  the  guasi-Lagrangian  coordinate  system  and  general- 
ized budget  equations  is  presented  in  Appendix  B.  The 
results  of  the  budget  analyses   are  presented  in  Chapter  IV. 


16 


Chapter   V    contains   the   conclusions     of   this    study   and    recom- 
mendations   for   further   research. 
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II.       DATA    BASE    AND    SYNOPTIC    OVERVIEW 

A.        DATA    BASE 

The  synoptic  descriptions  and  budget  analyses  of  the 
Presidents1  Day  Cyclone  are  based  solely  on  the  FG3E  Level 
Ill-b  data  set.  More  specifically,  FG3E  analyses  between 
1200  GMT  18  February  and  1200  GMT  20  February  1979  are  used 
to  define  the  development  and  movement  of  this  storm.  The 
FGGE  analyses  were  prepared  by  the  European  Center  for 
Medium  Range  Weather  Forecasts  (ECMWF)  using  a  three-dimen- 
sional optimal  interpolation  scheme  on  a  high  resolution 
(1.875  degrees  of  latitude  and  longitude)  nonstaggered  grid. 
This  nonlinear  analysis  scheme  produced  high-quality,  verti- 
cally and  horizontally  consistent  mass  and  momentum  fields. 
A  more  detailed  description  of  the  analysis  scheme  is  found 
in  the  Global  Weather  Experiment .  The  Daily,  Global  Analysis, 
£i£l  Jl,     (1981),    published   by   ECMWF. 

The  validity  of  the  FGGE  data  set  in  defining  the  Presi- 
dents* Day  Cyclone  was  demonstrated  by  Roman  (1981) .  His 
subjective  comparisons  with  the  Bosart  (1981)  analyses,  the 
Fleet  Numerical   Oceanography   Center  analyses   and  the   NMC   LFM 
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analyses    clearly   illustrate  the    accuracy    Df   the   FGGE   data    in 
representing   the   atmospheric   conditions   daring   the   period. 

Some  of  the  advantages  in  using  the  FGGE  data  set  as 
opposed   to   conventional   analyses    are: 

•  The  incorporation  of  a  variety  of  data  sources  includ- 
ing   off-time   and   late   data; 

•  The  high  resolution  regular  grid  is  particularly  useful 
in  examining  the  subsynoptic  scale  features  of  the 
Presidents'    Day   cyclone; 

•  The  data  set  includes  several  derived  fields  such  as 
relative   humidity   and   vertical  velocity; 

•  The  analyses  fully  encompass  the  area  of  development 
and    movement   of   the  Presidents'    Day   Cyclone. 

The  FGGE  Level  Ill-b  data  set  includes  fields  of  pres- 
sure height,  u  and  v  wind  components,  temperature  and  verti- 
cal velocity  for  all  standard  pressure  levels  and  relative 
humidity  for  all  standard  levels  up  to  300mb  (Table  I). 
Additionally,  absolute  vorticity  and  thickness  fields  were 
derived   for  use    in   this   thesis. 

B.       SYNOPTIC   OVERVIEW 

The  purpose  of  this  section  is  to  review  the  changes  in 
the  mass  and  momentum  fields  over  the  eastern  U.S.  and  west- 
ern Atlantic  that  led  to  the  development  of  the  Presidents' 
Day  Cyclone.  This  review  differs  from  previous  discussions 
of   the  storm    (Bosart,    1981;      and    Uccellini,      Kocin   and    Wash, 
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1981)  by  focusing  on  these  changes  in  relation  to  the  budget 
volume  translating  with  the  surface  pressure  minimum.  The 
discussion  is  based  on  the  pressure,  temperature  and  wind 
analyses  and  the  derived  vortic:ity  and  thickness  fields  from 
the  FGGE  data  set.  The  discussion  in  terms  of  vorticity  and 
thermal  advection  serves  to  provide  a  background  for  the 
circulation  and  angular  momentum  budget  study  analyses  pre- 
sented in  Chapter  IV. 

All  of  the  synoptic  figures  in  this  thesis  depict  budget 
volumes  with  lateral  boundaries  at  six  and  nine  degrees 
latitude.  The  six-degree  boundary  was  selected  because  it 
generally  includes  the  intense  part  of  the  storm.  The 
nine-degree  boundary  was  chosen  since  it  encompasses  nearly 
all  of  the  peripheral  features  related  to  the  storm.  The 
six  and  nine  degree  circles  are  referred  to  throughout  this 
section  as  the  inner  and  outer  boundaries.  The  circular 
representation  of  the  budget  volume  lateral  boundaries 
introduces  a  small  distortion  since  these  boundaries  should 
be  slightly  elliptical  on  the  conformal  projections  used  in 
the  figures.  Also,  the  position  of  the  budget  volume  at 
1200  GMT  18  February  is  based  on  the  first  distinguishable 
surface  pressure  minimum  of  the  Presidents1  Day  Cyclone  (not 
shown) . 
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TABLE    I 
FGGE    LEVEL    Ill-b    DATA    FIELDS    FOR    TIMES    SHOWN    BELOW 

SURFACE    PRESSURE    FIELDS     ARE    ALSO    AVAILABLE 


LEVEL    (mb    x    10)  =  100      85      70      50      40       30       25      20       15      10 
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U    WIND    COMPONENTS 
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X 

V    WIND    COMPONENTS 

X 

X 
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X 

X 

X 

X 

X 

X 

X 

TIME    PERIODS    AVAILABLE    (GMT) 


12  18  Jo  $6  12  18  dO  06  12 

18    FEB  19    FEB  20    FEB 


Figure  1  depicts  the  6-hourly  positions  of  the  storm's 
surface  pressure  minima  prior  to,  during  and  following  the 
period  of  explosive  cyclogenesis  from  3000  GMT  19  February 
through  0000  GMT  20  February  1979.  During  this  time,  the 
surface  pressure  minima  decreased  by  27  mb  for  an  average 
deepening   rate   in   excess   of   1    mb    per    hour   for    24   hours.      The 
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maximum  6-hour  decrease  was  10.3  mb  between  1200  GMT  and 
1800  GMT  19  February-  Figs.  2  through  6  depict  the  surface 
pressure  patterns  and  the  850  and  500mb  height  contours  and 
isotherms  at  12-hour  intervals  from  1200  GMT  18  February 
through    1200     GMT   20   February.  The    1000mb      thermal    fields 

appear  with  the  surface  pressure  patterns  since  the  surface 
temperature  fields  are  not  included  in  the  FGGE  Level-IIIb 
data   set.  The   500mb    wind     vectors   and      absolute   vorticity 

fields  at  12-hour  intervals  are  presented  in  Figs.  7  and  8. 
The  examination  of  these  figures  yields  qualitative  esti- 
mates 3f  mid-level  vorticity  advection.  Finally,  the 
1000-500mb  thickness  patterns  and  1000mb  height  contours 
appear  in  Figs.  9  and  10.  These  figures  are  useful  in 
determining   lower  tropospheric  thermal  advection. 

At  1200  GMT  18  February  the  dominant  surface  synoptic 
feature  is  a  massive  cold  high  pressure  system  centered  over 
northern  New  York  with  a  ridge  to  the  south  over  and  east  of 
the  Appalachians  (Fig.  2a)  .  This  ridge  is  a  hydrostatic 
reflection  of  the  shallow  cold  air  entrenched  east  of  the 
Appalachians.  The  thermal  trough  at  the  1000mb  level  in  this 
region  is  located  within  the  northwestern  portion  of  the 
outer  budget   volume      boundary.       Another      feature   of   interest 
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Figure  1.      Presidents      Day   Cyclone  track   at    6-hour   intervals 
from    1800   GMT   18   February   through    1200   GMT   20 
February,    1979   with   5    degree   radius    circles. 
Times   are   date   and   hour   in   GMT    (Z)  .    Central 
surface   pressures   in    mb  appear   below   times. 
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at  the  surface  is  the  inverted  pressure  trough  in  the 
central  Gulf  of  Mexico.  Additionally,  a  *eak  trough  is  seen 
near  the  center  of  the  budget  volume.  The  850  mb  contour 
pattern  (Fig.  2b)  shows  a  strong  ridge  over  the  northeastern 
states  that  extends  into  the  Great  Lakes  region,  and  a  sharp 
trough  from  North  Dakota  southeastward  across  the  inner 
budget  volume  boundary.  The  85  0mb  isothsrm  pattern  is  pri- 
marily zonal  in  the  vicinity  of  ths  budget  volume  with  lit- 
tle indication  of  the  underlying  shallow  cold  air.  The  key 
feature  at  the  500mb  level  (Fig.  2c)  is  a  short  wavs  embed- 
ded in  a  nearly  zonal  flow  across  the  upper  Midwest.  The 
wave  is  generally  in  phase  with  the  thermal  trough  at  this 
level.  A  weak  ridge  is  seen  in  the  503mb  contour  pattern 
north  of  the  inner  volume  boundary.  The  absolute  vorticity 
maximum  (Fig.  7a)  extends  into  the  northern  portion  of  the 
budget  volume  and  southward  into  Texas.  Positive  vorticity 
advection  (PVA)  is  occurring  to  the  west  and  northwest  of 
the  budget  volume  from  Illinois  southwest  ward  into  Louisi- 
ana. The  1000-500mb  thickness  pattern  (Fig.  9a)  indicates  a 
broad  area  of  strong  cold  advaction  throughout  the  lower 
troposphere  over  the  western  North  Atlantic.  Weak  cold 
advection  is  present  across  the  northeastern  budget  volume 
boundary. 
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Figure  2.      a.    Surface   pressure  and  10Q0mb   isotherms,    b. 

850mb   height   contours   and    isotherms   and   c.    SOOmb 
contours   and    isotherms   for    1200    GMT    18    February 
1979.    Surface   pressure   is   p-1000rab.      Contours   are 
labeled  in   decameters    and    isotherms    (°C)    are 
dashed.      Circles  are    6   and    9   degree   budget    volume 
lateral   boundaries. 
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Significant  changes  occur  in  the  low-level  patterns 
during  the  12  hour  period  ending  at  0000  GMT  19  February. 
The  surface  high  pressure  center  (Fig.  3a)  moves  southeast- 
ward to  Massachusetts.  More  importantly,  however,  is  the 
development  of  a  cyclonic  circulation  off  the  coast  of  South 
Carolina.  There  was  little  indication  12  hours  earlier  that 
cyclogenesis  would  commence  in  this  area.  One  might  argue 
that  the  inverted  trough  previously  d,n  the  Gulf  could  have 
signaled  the  onset  of  cyclogenesis,  if  the  trough  had  moved 
to  the  current  position  of  the  low  with  the  500mb  flow. 
However,  the  1800  GMT  surface  analysis  (not  shown)  shows  no 
indication  of  the  trough  translating  northeastward  across 
the  Florida  peninsula.  Coinciding  with  the  cyclonic  devel- 
opment off  the  coast  is  a  reorientation  and  strengthening  of 
the  1000mb  thermal  gradient  along  the  coast  from  Georgia  to 
Maryland  within  the  inner  volume  boundary.  The  presence  of 
the  coastal  front  and  its  contribution  to  the  cyclone  devel- 
opment is  discussed  by  Bosart  (1981).  In  addition,  a  sig- 
nificant thermal  ridge  has  developed  within  the  inner  volume 
boundary. 

The  primary  change  in  the  850mb  contour  pattern  (Fig. 
3b)   is  the   southeastward  extension  of  the   trough  over  the 
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surface  low.  This  trough  now  sr.compasses  the  entir=  north- 
west quadrant  of  the  budget  volume.  Little  change  is  seen 
in  the  850mb  isotherm  pattern  during  the  12-hour  period. 
The  thermal  patterns  at  the  130  0  and  353mb  levels  indicate 
low-level  destabilization  within  the  innsr  volume  boundary 
which  would   favor  organized  convection. 

The  500mb  wave  (Fig.  3c)  has  moved  eastward  to  Illinois 
and  is  more  in  phase  with  the  thermal  trough  at  this  time. 
The  wave  is  beginning  to  penetrate  the  northwest  quadrant  of 
the  cuter  budget  volume  boundary.  A  weak  ridge  is  seen  in 
the  isotherm  pattern  at  this  lsvel  throughout  the  inner  vol- 
ume boundary.  The  vorticity  pattern  (Fig.  7b)  shows  weak 
PVA  over  the  northern  half  of  the  volums.  Thermal  advec- 
tion, as" depicted  in  the  thickness  pattarn  (Fig.  9b),  is 
confined  to  the  inner  volume  radius  with  warm  advection  to 
the  east  and  weak  cold  advection  to  the  wast  of  the  cyclone 
center. 

The  1200  GMT  19  February  surface  pressure  pattern  (Fig. 
Ua)  indicates  significant  cyclo genesis  has  occurred.  The 
minimum  pressure  decreased  by  9  mb  in  the  12-hour  period 
ending  at  this  time.  The  surface  isobars  show  a  highly  sym- 
metrical  circular   cyclone   structure      located   entirely    within 
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Figure  3.      Same   as   Fig.    2   except    for    0000    GMT    19   February 
1979. 
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the  inner  budget  radius.  The  100Dmb  thsrmal  ridge  is  con- 
siderably sharper  and  is  in  nearly  perfact  alignment  with 
the  cyclone  center  and  associated  pressure  trough.  The 
1000mb  isotherm  pattern  reveals  the  devslopment  of  a  baro- 
clinic  zone  to  the  southwest  of  the  cyclone  within  the  inner 
budget  volume.  The  850mb  height  contours  (Fig.  4b)  now  show 
closed  cyclonic  circulation  centered  in  the  northwest  quad- 
rant of  the  inner  volume  radius  which  indicates  a  northwest 
tilt  of  the  cyclone  axis  toward  the  cold  air.  Significant 
changes  in  the  850mb  isotherm  pattern  have  occurred  during 
this  12-hour  period.  The  development  of  a  thermal  ridge 
over  the  surface  center  indicates  that  ths  warm  core  is  pen- 
etrating to  higher  levels.  This  warming  is  liksly  the 
result  of  latent  heat  release  through  oonvection  and  the 
warm  advection  shown  at  the  850m b  level  12  hours  earlier. 

The  500mb  wave  (Fig.  4c)  is  now  centered  over  Pennsylva- 
nia and  is  nearly  in  phase  with  the  thermal  trough  within 
the  northwest  quadrant  of  the  outer  budget  volume  boundary. 
The  500mb  vorticity  center  (Fig.  7c)  is  now  within  the  outer 
radius  with  significant  PVA  in  the  inner  radius  primarily  to 
the  northeast  of  the  surface  center.  The  thickness  pattern 
(Fig.   9c)   shows  a  warm  ridge   centered  on  the  cyclone  with 
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warm  advection  throughout  the  eastern  half  and  cold  advec- 
tion  throughout  the  western  half  of  the  outer  budget  radius. 
Intense  cyclogenesis  continues  during  the  next  12-hour 
period  ending  at  0000  GMT  20  ?abruary.  During  this  period 
the  surface  pressure  minimum  (Pig-  5a)  decreases  by  nearly 
18  mb.  At  the  same  time,  the  cyclone  takes  a  more  easterly 
course  following  the  maximum  500mb  PVA  region.  The  vertical 
structure  indicates  that  the  cyclone  is  nearing  the  mature 
stage.  The  circular  surface  and  950mb  centers  are  nearly 
coincident  and  the  isotherm  patterns  show  a  vertical  align- 
ment of  the  thermal  ridge  through  the  850ab  level  (Fig.  5b). 
The  500mb  wave  (Fig.  5c)  is  almost  directly  over  the  surface 
cyclone  and  is  encompassed  by  the  outer  budget  radius.  The 
thermal  trough  now  lags  the  500mb  wave  with  an  intensified 
thermal  ridge  within  the  eastern  half  of  the  outer  volume 
boundary.  Thus,  the  warm  cora  now  extends  through  the  low 
troposphere.  The  thickness  pattern  (Fig.  10a)  clearly 
illustrates  a  mature  warm  core  structure.  The  vertical 
extension  of  the  warm  core  through  the  500mb  level  is  most 
likely  due  to  convection  since  little  thermal  advection  is 
seen  at  this  level  prior  to  0000  GMT.  The  mature  structure 
indicates  that  little   further  deepening  of  the   cyclone  can 
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Figure  4.      Same   as   Fig.    2   except    for    1200    GMT    19    February 
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be  expected.  The  500mb  vorticity  pattern  (Fig-  8a)  shows 
PVA  throughout  the  eastern  half  of  the  budget  volume  with 
maximum  PVA  due  east  of  the  cyclone  center. 

The  surface  pressure  pattern  (Fig.  5a)  at  1200  GMT  20 
February  shows  that  the  cyclone  center  has  accelerated  in  an 
easterly  direction  during  the  previous  12  hours  with  little 
further  deepening.  The  storm  center  moves  nearly  250  nm 
eastward  for  an  average  velocity  of  approximately  20  knots. 
Little  change  is  seen  in  the  thermal  structure  through  the 
low  troposphere  in  the  inner  budget  volume.  This  might 
indicate  that  a  moist  adiabatic  lapse  rate  has  been  achieved 
by  convection  throughout  much  of  the  inner  budget  radius. 
The  500mb  vorticity  pattern  (Fig.  9b)  now  shows  moderate  to 
strong  PVA  in  the  eastern  half  and  definite  negative  vortic- 
ity advection  (NVA)  in  the  northwest  quadrant  of  the  outer 
volume  boundary. 

The  intent  of  this  synoptic  overview  has  been  to  present 
the  highlights  of  the  development  of  the  Presidents'  Day 
Cyclone  in  terms  of  the  mass  and  circulation  fields  as  a 
prelude  to  the  storm-following  budget  analyses  in  this  the- 
sis. The  goal  of  budget  studies  of  this  storm  is  to  quan- 
tify  the    forcing   mechanisms   of   explosive   maritime 
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cyclogenesis.  A  brief  discussion  of  the  quasi-Lagrangian 
technique  employed  in  the  budget  study  is  presented  in  the 
next   chapter. 
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Figure  5.   fame  as  Fig.  2  except  for  0000  GMT  20  February 
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Figure  6.      Same   as   Fig.    2   except    for    1200    GMT   20   February 
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Figure  7.      500mb  absolute    vorticity   fields    (x    10-s   sec~M 
^S3oW12d   ?fS£ors  for    a-    1200   SMT    18   February 
l219*nbZ   2000   GMT   19    February   1979   and  c.    1200 
GMT   \9   February    1979.    Circles   are   6    and   9    degree 
budget   volume   lateral    boundaries. 
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Figure  8.      Same   as   Fig.    7    except    for    a.    0300    GHT   20   February 
*  1979    and   b.    1200  GMT    20   February    1979. 
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Figure  9.  1000-500mb  thickness  and  lOOQmb  contour  patterns 
for  a.  1200  GMT  T8  February  1579.  b.  0000  GMT  19 
February  1979  and  c.  1200  GMT  19  February  1979. 
Thickness  contours  are  dashed  and  labeled  in 
decameters-  1000  mb  contours  are  labeled  in 
meters.  Circles  are  5  and  9  degree  budaet  volume 
lateral   boundaries. 
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Figure  10.   Same  as  Fig.  9  except  a.  0000  GMT  20  February 
1979  and  b.  1200  GMT  20  February  1979.       J 
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III.       QUASIzLAGRANGIAN    DIAGNOSTIC    TECHNIQUE 

The  objective  of  budget  studies  is  to  evaluate  the 
changes  of  atmospheric  properties  within  a  specific  atmos- 
pheric volume-  Normally  these  studies  have  been  used  prima- 
rily to  study  the  general  circulation  and  planetary  scale 
balances  of  atmospheric  properties  such  as  mass,  momentum 
and  kinetic  energy.  For  synoptic-scale  studies  the  approach 
has  been  to  evaluate  changes  in  these  properties  within  a 
stationary  volume  (Rung  and  Baksr,  1975;  Smith,  1973).  The 
effects  of  synoptic-scale  disturbances  migrating  through  the 
volume  are  quantitatively  evaluated  relative  to  the  fixed 
volume  only.  This  approach  is  not  well  suited  to  the  study 
of  cyclone  development  since  changes  in  the  budget  volume 
can  be  due  to  the  movement  of  the  cyclonic  circulation  into 
or    out  of   the  volume   as    well  as    due   to   cyclone   development. 

The  QLD  budget  approach  employed  in  this  thesis  uses  a 
budget      volume   which      translates      with     the   cyclone.  This 

method  permits  a  quantitative  evaluation  of  the  changes  in 
the  budget  property  relative  to  the  cyclone.  Specifically, 
this  system  isolates  the  horizontal  advection  and  transport 
relative  to  the  cyclone  from  the  horizontal  advection  and 
transport    associated      with    the      translation    of      the  cyclone. 
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Further,  the  vertical  redistribution  of  the  budget  property 
is  the  result  of  the  vertical  transport  processes  in  the 
vicinity  of  the  cyclone.  Thus,  the  QLD  technique  is  well 
suited  to  the  study  of  cyclogenesis.  Table  II  presents  an 
updated  list  of  recently  completed  studies  employing  the  QLD 
method. 


TABLE  II 
PREVIOUS  QUASI-LAGRANGIAN  DIAGNOSTIC  BUDGET  STUDIES 

(EXTRACTED  IN  PART  FROM  WASH,  1978) 

Property  Researcher  (s) 

Available  Potential  Energy      Spaete  1974 
Circulation  Wash  1975 

Mass  Johnson  and  Downey  1976 

Absolute  Angular  Momentum       Johnson  and  Downey  1976 
Kinetic  Energy  Chen  and  Bosart  1977 

Absolute  Angular  Momentum       Wash  1978 
Mass  Roman  1981 

Mass  Tallman  1982 


The  budget  equation  is  evaluated  with  respect  to  a 
budget  volume  in  spherical  cosrdinates  (Appendix  B) .  In 
this  thesis,  the  budget  volume  is  centered  over  and  main- 
tains its  position  relative  to   the  surface  pressure  minimum 
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of  the  Presidents'  Day  Cyclone.  The  volume  is  defined  in 
the  vertical  by  10  mandatory  pressure  lsvels  (1000,  850, 
700,  500,  400,  300,  250,  200,  150  and  100  mbs)  and  horizon- 
tally by  radii  for  example,  six  degrees  of  latitude.  Appen- 
dix B  establishes  the  quasi-Lagrangian  budget  volume 
coordinate  system  and  the  generalized  budget  equations  that 
can    be  applied  to   any   selected   atmospheric    property. 

The  first  step  in  the  calculation  of  the  budgets  is  to 
interpolate  the  u  and  v  wind  component  values  at  the  rectan- 
gular gridpoints  in  the  FGGE  data  set  to  the  budget  volume 
gridpoints.  The  interpolated  winds  are  then  converted  to 
normal  and  tangential  components  relative  to  the  circular 
budget  volume  boundaries.  Additionally,  the  movement  of  the 
budget  volume  (W0)  is  used  to  compute  the  normal  and  tangen- 
tial  components   relative  to  the    translating    volume. 

The  budget  equations  for  the  circulation  and  angular 
momentum  budgets  in  pressure  coordinates  are  presented  in 
Tables  III  and  IV.  In  the  vorticity  budgst  formulation,  the 
horizontal  advection  and  divergence  terms  are  combined  to 
form      the      lateral      transport      using      the      Gauss      Divergence 

Theorem,  f   /-(u-w)-vc   -?   v«(u-w)dA 

~*  ~  a      a         "w  ~ 

A 

—  /    /     V*C     (U-W)dA    =-/C     (U-W).dL. 
a    "■*  ~~  a    ~  ~    3 

A 
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The  vertical  redistribution  occurs  through  the  vertical 
advection  term  and  the  friction  and  tilting  terms  complete 
the  budget.  The  angular  momentum  budget  is  formulated  by 
combining  the  continuity  and  angular  momentum  equations. 
Since  there  is  no  divergence  term,  the  vertical  redistribu- 
tion occurs  through  a  vertical  transport  term.  The  friction 
term  completes  this  budget. 

The  budget  volume  grid  consists  of  36  equally  spaced 
points  around  each  circular  lateral  boundary  at  each  pres- 
sure level.  Absolute  vorticity  and  angular  momentum  values 
are  computed  at  each  of  the  volume  gridpoints.  The  calcula- 
tions of  the  time  tendency,  transport  and  source/sink  terms 
are  based  on  these  gridpoint  values.  The  tilting  term  in 
the  vorticity  equation  is  computed  as  a  source.  There  are 
no  sources  of  absolute  angular  momentum  within  the  budget 
volume,  and  the  only  sink  is  due  to  friction  which  is  com- 
puted using  a  stability  independent  parameterization  scheme 
(Johnson  and  Downey,  1976)  with  a  constant  geostrophic  drag 
coefficient  (0.035). 

The  transport  terms  are  broken  into  msan  and  eddy  modes 
(Table  7).   The   purpose   of  this  statistical  partitioning 
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is  to  isolate  changes  in  the  budget  property  associ- 
ated with  the  synoptic-scale  mean  and  sidy  circulations. 
The  mean  and  eddy  modes  represent  the  transport  by  the  irro- 
tational  and  rotational  wind  components,  respectively.  In 
this  context,  the  mean  and  eddy  modes  ace  analogous  to  the 
divergence  and  horizontal  advection  terms  in  the  vorticity 
eguation.  In  this  study  the  divergence  and  advection  terms 
are  calculated  for  comparison  with  the  mean  and  eddy  modes. 
The  comparisons  serve  to  further  substantiate  the  validity 
of  the  mean  and  eddy  mode  partitioning  in  the  QLD  technique. 
The  FGGE  vertical  velocity  fields  were  calculated  based 
on  divergences  derived  through  a  nonlinear  normal  mode  ini- 
tialization scheme.  An  objective  of  this  thesis  is  to  com- 
pare the  FGGE  vertical  velocities  with  ths  vertical  velocity 
values  derived  via  the  QLD  method  in  the  mass  budget  of  the 
Presidents1  Day  Cyclone  (Roman,  1981) .  The  purpose  of  this 
comparison  is  to  aid  in  the  assessment  of  the  vertical 
velocities  computed  in  the  budget  technique.  The  comparison 
is  presented  in  conjunction  with  the  analysis  of  the  budget 
vertical  transport  term  in  the  next  chapter. 
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TABLE    III 

QUASI-LAGRANGIAN    CIRCULATION    BUDGET    EQUATION    IN    ISOBARIC 

VERTICAL    COORDINATES 

(after  Wash,    1978) 


lOOmb      8       2ir    , 
C   -      f  f      f         -C    r2sin6    dadB(-dp) 

a    lOOOmb    0       0  9    a 


where  C     is   absolute   circulation  and  £     is 
a  *a 

absolute  vorticity. 


The  budget  equation  is 


|— a=LT(C    )+VT(C    )+S(C    ) 
6 1  a  a  a 


where  the  lateral  transport  is 

lOOmb      2it    -, 
LTtC    )=      /  /         -(U-W)ftC    rsinS    da(-dp) L 

a      lOOOmb    0         ^  »a  0 


and  the  vertical  advection  is 


VTCC    )=    /      /         "d>:^a    r2sinS    dad6. 
a         0      0         9    3P 


The  source/sink  term  is 

lOOmb      8       2tt  1    dc         9  „  ,     .    * 

SCC    )-      /  /       /  -  5£a    r2sin8    dodB(-dp), 

a       lOOOmb    0       0  g    at 


See  Appendix  B  for  notation  explanation. 
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TABLE    IV 

OUASI-LAGEANGIAN    ANGULAR    MOMENTUM    BUDGET    EQUATION    IN 
y  ISOBARIC    VERTICAL    COORDINATES 

(after  Mash,    1978) 


lOOmb      8       2ff    i         . 
6  -      /  /      /        -g    r2sinB    dadB(-dp) 

a    lOOOmb    0      0         g    a 


where  G     is  the  volume  integral  of  absolute  angular 
momentum  g  . 

SI 


The  budget  equation  is 

«G 

-r^a-LTtG    )+VT(G  J  +  S(G    ) 

o  t  a  a  a 


where  the  lateral  transport  is 
LT 


lOOmb       2tt    I  . 

(G   )-      /  /        -<£-&)    grsine    do(-dp)| 

a      lOOOmb    0         9  p 


and  the  vertical  divergence  of  the  vertical  transport  is 

lOOmb      3      2ir    la 
VT(G   )«      /  /      /        ~-(wg    )r2ain3    dodB  (-dp). 

a      lOOOmb    0       0         *3P         a 

The  sink  term  is 

lOOmb      S       2»    i    da         , 
S(g    )-      /  /      /        *■  TT»    r2sinB    dadS(-dp). 

a       lOOOmb    0       0         9    at 


See  Appendix  B  for  notation  explanation. 
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TABLE  V 
PARTITION  OF  LATERAL  TRANSPORT  INTO  MEAN  AND  EDDY  MODES 

(after  Wash,  1978) 


lOOmb   2ir  i  _ ;  i 

EM(C  )=   /     /    -  C*(U-W)*  rsinB  da(-dp)|a 
a   lOOOmb  0    ^   a  -  -  6  9 


lOOmb   2ir  i    _a 


MM(C  )•   /     /    -  C   (U-W)ft  rsinB  da(-dp)|ft 
a   lOOOmb  0    9   a       p  p 


lOOmb       2tt    i    ,  , 

M(G    )=       /  /         -   gM£-W)*    rsinB    da(-dp) 


lOOmb       2ir    i 
OOOmb    0         9 

lOOmb     ^2ir    j_    _a    a 


MM(G    )=      /  /         -  gZ    CU-W)£    rsinB    do(-dp)L 

a       lOOOmb    0         9       a  s  * 


EH  is  the  eddy  mode   and' MM  is  the  mean  mode. 

-   a       1         2v 
Note  that         (    )    =   :r-  /         (    )da 

2tt    Q 


and         ()**()-()  °      is  the  deviation  of 


the  property   from  its  mean  around  a  lateral 
boundary.      See   Tables   III   and  IV  and  Appendix 
B   for  notation  explanation. 
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IV.       SDASIzLAGRANGIAN    BUDGET    ANALYSIS 

A.        GENERAL 

The  absolute  circulation  or  vorticity  and  absolute  angu- 
lar momentum  budget  results  are  presented  in  this  chapter. 
Emphasis  is  placed  on  the  vorticity  budget  since  the  results 
are  similar  and  vorticity  analyses  are  mors  commonly  used  to 
describe  cyclone  development.  The  angular  momentum  results 
are  discussed  in  depth  where  significant  differences  exist, 
particularly   in   the   vertical  transport   term. 

The  results  are  presented  by  means  of  layer  or  level/ 
time  or  time  period  cross-sections  and  selected  vertical 
profiles.  Results  are  presented  for  two  radii,  6°  and  9°, 
which  were  illustrated  in  the  synoptic  discussion.  Values 
are  displayed  on  levels  or  layers  over  6-hour  time  periods. 
For  example,  the  600mb  layer  extends  frDm  700  to  500mb  and 
the  1500  GMT  time  period  runs  from  1200  to  1800  GMT.  The 
sign  convention  is  positive  for  processes  which  produce  a 
vorticity  or  angular  momentum  increase  such  as  inward  trans- 
port, inward  horizontal  advection  and  sources.  Upward  ver- 
tical velocity    (omega)    is  negative    (pressure    coordinates). 
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B.   TIME  TENDENCY 

The  dramatic  buildup  of  absolute  vorticity  within  the 
budget  volume  during  the  development  of  the  Presidents'  Day 
Cyclone  is  illustrated  by  the  vorticity  time  sections  in 
Fig.  11.  These  patterns  rsflecx  the  vorticity  maximum 
developing  within  the  budget  radii  at  the  1000  mb  level 
(Fig.  12)  and  the  vorticity  maxima  moving  into  the  budget 
circles  at  the  500  and  300  mb  levels  (Figs.  7  and  13 
respectively) . 

The  vorticity  time  section  patterns  (Fig.  11)  within 
both  radii  are  similar  with  minima  occurring  in  the  225mb 
layer  during  the  1500  GMT  18  February  time  period.  These 
minima  are  seen  in  terms  of  anticyclonic  shear  and  carvature 
within  the  budget  radii  in  the  3  00mb  wind  field  (Fig.  14a) . 
The  vorticity  values  are  greatar  in  the  6°  radius  than  those 
in  the  9°  radius  after  the  0900  GMT  19  February  period.  The 
larger  gradient  with  respect  to  time  in  the  inner  radius 
depicts  the  rapid  increase  and  concentration  of  cyclonic 
circulation  as  the  storm  develops.  Maxima  in  both  volumes 
occur  in  the  upper  layers  during  the  0930  GMT  20  February 
period.  The  larger  vorticity  increase  with  time  in  the 
upper  layers  is  due  to  the  development  of  cyclonic  shear  and 
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curvature  as  the  upper  tropospheric  trough  and  associated 
jet  maximum  dominate  the  southern  half  of  the  inner  budget 
volume  (Fig.  14b).  The  vorticity  profiles  (Fig.  15)  at  the 
beginning  and  end  of  the  intense  cyclogenesis  period  show 
that   the    largest   increase   occurs    in   the   upper    layers. 

Vorticity  time  tendencies  computed  over  6-hour  intervals 
are  depicted  in  Fig.  16.  These  patterns  show  positive  ten- 
dencies in  all  layers  and  time  periods  with  the  exception  of 
the  lower  layers  within  the  inner  radius  during  the  final 
time  period.  Both  patterns  show  primary  and  secondary  max- 
ima. The  primary  maximum  in  the  9°  radius  (350mb  layer, 
0300  GMT  19  February  period)  occurs  earlier  than  the  corre- 
sponding maximum  in  the  6°  radius.  The  outer  region  of  the 
cyclone  is  spinning-up  earlier  in  response  to  the  influence 
of  the  upper  level  trough  and  jet  maximum  at  this  time.  As 
the  trough  overtakes  the  lower  level  cyclone,  the  6°  radius 
shows     a   large      vorticity   increase.  The    secondary      maxima 

occur  during  the  same  period  (2100  GMT  19  February)  in  each 
volume,  but  at  a  higher  layer.  These  maxima  are  associated 
with  the  strong  cyclonic  shear  and  curvature  north  of  the 
upper  level  jet  core  which  is  located  within  the  budget  vol- 
ume   at  this  time    (Fig.    14b).      The   magnitude    of   the   secondary 
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Figure  11. 


Absolute    vorticity   time  sections       (x   10~5 
sec-M- 
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Figure  12. 


1000mb  wind   vectors    and  absolute    vorticity    (x 
10-s   sec-M     for  a.    3  000   GMT   19    February   1979. 
1200   GMT   19    February    1979   and   c.      0000   GMT    20 
February    1979.      Circles  are  6    and   9    degree 
budget   volume   lateral   boundaries. 
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Figure   13.      300    cab   wind   vectors 

(x   10-s   sec-M    for   a.    0600   GMT 

b.    1800   GMT    19    February   1979.      Circles    are   6    and 

9   degree    budget  volume  lateral    boundaries. 


and   and   absolute   vorticity 
.    0600   GMT    19   February    and 
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Figure   14.      300    mb   wind   vectors    and   isotachs    (m  sec-M     for 
a.    1800   GMT    18   February   an8  5.0600   GMT  20 
February    1979.      Circles    are   6    and   9    degree 
budget   volume   lateral   boundaries. 
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Figure  15. 


125  — 
175  |— 

225 


19002     2000Z 


JJ350 

5 

W 

a  45° 

600 
775 
925 


5  10         15 
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Absolute  vorticity  profiles  for  radius  six 
(solid)  and  radius  nine  (dashed) . 


maximum  in  the  6°  radius  is  nearly  twice  that  in  the  9° 
radius,  because  the  maximum  cyclonic  curvature  and  shear 
vorticity  occur  within  the  inner  boundary. 

The  absolute  angular  momentum  time  section  patterns 
(Fig.  17)  are  similar  to  those  of  absoluts  vorticity,  except 
the  values  in  the  outer  radius  exceed  those  in  the  inner 
radius.  This  difference  is  because  storn  angular  momentum, 
ga  =Jc  •  R  x  (&-WL,  is  proportional  to  the  distance  (R)  from 
the  storm  center.    The  absolute  angular  momentum  (Fig.   18) 
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Figure   16. 


Absolute    vorticity   time   tendencies   over   6    hours 
(x    10-l°    sec-2).      Dashed  indicates   negative 
values. 
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and  vorticity  time  tendency  patterns  exhibit  remarkable  sim- 
ilarities with  virtually  coincident  maxima.  The  larger 
angular  momentum  tendencies  in  the  larger  volume  again 
reflect  the  dependence  of  this  property  on  radius.  The 
buildup  of  angular  momentum  within  the  budget  volume  por- 
trayed in  these  patterns  is  consistent  with  the  findings  of 
Johnson  and  Downey  (1976). 

In  summary,  the  vorticity  and  angular  momentum  patterns 
show  that  the  largest  increases  occur  within  the  24-hour 
period  of  intense  cyclogenesis  (0000  GMT  19  February  through 
0000  GMT  20  February).  The  maximum  vorticity  buildup  occurs 
in  the  upper  layers  because  the  polar  jet  core  moves  into 
the  budget  volume.  The  primary  difference  between  the  vor- 
ticity and  angular  momentum  buildup  is  that  the  maximum 
angular  momentum  is  found  at  a  lower  layer  particularly  in 
the  6°  radius  volume  (Fig.  17). 

C.   LATERAL  TRANSPORT 

The  lateral  transport  of  absolute  vorticity  (Fig.  19)  is 
the  most  important  forcing  term  in  the  budget.  The  maximum 
inward  transport  occurs  in  the  upper  layers,  although  the 
transport  is  positive  for  nearly  all  of  the  layers.  The 
maximum  in  the   smaller  volume  occurs  later,   is   in  a  lower 
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layer  and  is  considerably  greater  in  magnitude  than  the  max- 
imum in  the  larger  volume.  Both  maxima  occur  within  the 
24-hour  period  of  intense  cyclDgenesis;  in  particular,  the 
larger  maximum  in  the  inner  volume  occurs  during  ths  period 
of  most  rapid  deepening  (1200-1800  GMT  19  February).  The 
maximum  transport  in  the  larger  volume  occurs  earlisr  (0600 
GMT  19  February)  as  the  vorticity  maxinum  approaches  the 
northwestern  portion  of  the  outer  boundary  (Fig.  13a) . 
Twelve  hours  later,  the  vorticity  maxiaum  crosses  the  6° 
boundary  (Fig.  13b)  and  produces  the  larger  maximum  trans- 
port at  the  inner  radius. 

As  described  in  the  preceding  chapter,  the  total  lateral 
transport  can  be  partitioned  into  mean  and  eddy  mode  compo- 
nents. The  purpose  of  this  partitioning  is  to  break  t-he 
total  transport  into  components  to  improve  the  interpreta- 
tion of  the  lateral  transport  dynamic  forcing.  The  mean 
mode  represents  the  mean  vorticity  transport  by  the  average 
normal  wind  component  (mean  divergence  or  convergence) , 
while  the  eddy  mode  is  the  transport  due  to  the  covariance 
of  the  vorticity  and  normal  wind  component  deviations.  In 
more  general  terms,  the  mean  mode  transport  is  due  to  the 
mean  inflow  and  outflow   of  the  cyclone  and  the   eddy  mode 
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Figure   18.      Absolute   angular  momentum  time   tendencies    (x    10s 
cm2    sec*2)  .      Dashed    indicates   negative   values. 
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Figure   19. 


Absoluts    vorticity   lateral  transport    (x   10~l° 
sec-2) .      Dashed  lines  indicate    negative   values 
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transport  reflects  transport  due  to  the  asymmetries  in  the 
cyclone  and  jet  streak. 

An  alternative  partitioning  of  the  lateral  transport 
term  follows  the  traditional  vorticity  equation  with  the 
horizontal  advection  and  divergence  terms.  The  area-aver- 
aged absolute  vorticity  horizontal  advection  was  calculated 
for  comparison  with  the  eddy  mode  transport.  The  divergence 
term  was  not  calculated  in  this  study  and  its  comparison 
with  the  mean  mode  transport  must  be  accomplished 
indirectly. 

The  vorticity  eddy  mode  transport  time  sections  (Fig. 
20)  are  very  similar  to  the  total  transport,  except  the  mag- 
nitudes are  larger.  The  maximum  in  the  smaller  volume  again 
occurs  during  a  later  time  period  and  is  greater  in  magni- 
tude than  the  outer  radius  maximum.  The  primary  features  in 
the  area-averaged  horizontal  vorticity  advection  patterns 
(Fig.  21)  generally  correspond  to  those  in  the  eddy  mode 
patterns.  The  upper  level  maxima  in  the  inner  volume  occur 
during  the  same  time  period,  although  the  PVA  maximum  occurs 
in  a  higher  layer.  The  large  maximum  in  the  upper  tropo- 
sphere results  from  the  movement  of  the  strong  vorticity 
maximum  and  jet  streak  into  the  cyclone  region  as  described 
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in  the  synoptic  discussion  (Figs.  7,  8,  13  and  14).  The 
similarities  in  these  patterns  suggest  a  strong  correlation 
between  the  eddy  mode  and  the  traditional  view  of  the  role 
of  PVA  at  upper  levels  in  the  development  of  extratropical 
cyclones    (Petterssen,    1956). 

The  most  significant  difference  between  the  eddy  mode 
and  area-averaged  advection  patterns,  particularly  within 
the  inner  volume,  is  the  PVA  (approximately  15  x  10-*° 
sec-2)  in  the  925  rab  layer  during  the  2100  GMT  19  February 
period.  This  low-level  PVA  does  not  agree  with  the  advec- 
tion implied  by  the  1000mb  wind  vectors  and  absolute  vortic- 
ity  patterns  (Fig.  12).  These  patterns  actually  indicate  a 
net  NVA  within  the  inner  radius  since  the  wind  vectors  are 
par'allel  or  cross  the  vorticity  isopleths  toward  higher  val- 
ues. In  an  effort  to  resolve  this  discrepancy,  the  absolute 
vorticity  advection  values  at  each  gridpoint  were  examined. 
The  data  indicate  that  the  low-level  PVA  is  occurring  prima- 
rily at  the  inner  radii  (less  than  6°)  to  the  east  and 
northeast    of     the  storm   center.  The   low-level   PVA      in   the 

time  sections  is  the  result  of  asymmetries  in  the  low-level 
vorticity  contours  and  winds  relative  to  the  translating 
storm   volume.        The     effect  of  the  translating      volume    is   to 
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enhance  the  low-level  advecting  (relative)  winds  (easter- 
lies) and  to  reduce  the  upper  level  advecting  winds  (wester- 
lies) .  This  accounts  for  the  smaller  difference  in  magni- 
tude between  the  upper  and  lower  level  vorticity  advection 
maxima  in  a  quasi-Lagrangian  compared  to  a  Eulerian  budget 
study-  Aside  from  these  small  discrepancies,  the  overall 
similarities  in  the  patterns  support  the  close  relation 
between  the  eddy  mode  lateral  transport  and  horizontal 
advection. 

The  mean  mode  vorticity  patterns  (Fig.  22)  show  outward 
lateral  transport  in  the  mid  to  upper  troposphere  and  inward 
transport  in  the  lower  layers.  These  patterns  reflect  vor- 
ticity transport  due  to  the  low -level  convergent  and  upper- 
level  divergent  cyclone  structure.  The  upper-level  outflow 
in  the  smaller  volume  increases  with  time  to  the  end  of  the 
period.  In  contrast,  the  smaller  maximum  outflow  value  for 
the  larger  volume  is  achieved  relatively  early  and  is  main- 
tained thereafter.  The  larger  mean  mode  inflow  within  the 
6°  radius  volume  is  the  result  of  larger  vorticity  values 
within  the  inner  boundary  as  well  as  the  more  vigorous  mass 
circulation   at  the  inner  radii. 


64 


18FEB 
\25r 


19FEB 


2DFEB 


9°    VOLUME 


6°    VOLUME 


03Z     09Z     152     21Z     03Z 
t:«e  PER I BOS 


09Z 


Figure  20. 


Absolute    vorticity   eddy   mode   lateral   transport 
(x    10-*o    sec-2). 
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Figure  21. 


Areal   average   absolute  vorticity   horizontal 
advection    (x    10-l<>   sec2)  . 
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Figure  22. 


Absolute    vorticity  i ean   mode   lateral  transport 
(x    10~10    sec-2).      Dashed    lines   indicate  negative 
values. 
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The  vorticity  lateral  transport  profiles  (Fig.  23) 
illustrate  the  mean  and  eddy  mode  contributions  to  the  total 
lateral  transport.  The  mean  mode  outflow  (corresponding  to 
divergence)  is  overwhelmed  by  the  eddy  mode  inflow  (PVA)  in 
the  upper  layers  while,  in  the  lower  layers  the  msan  mode 
inflow    (convergence)       dominates.  The   relationship   between 

the  mean  mode  transport  and  the  divergence  is  also  illus- 
trated by  noting  that  the  level  of  the  zero  value  in  the 
mean  mode  profile  (690mb)  agrees  with  the  mid-tropospheric 
level  of  nondivergence  (LND)  implied  by  the  lateral  mass 
transport   patterns    (Roman,    1981)  . 

The  absolute  angular  momentum  lateral  transport  patterns 
(Figs.  24-26)  are  also  similar  to  the  corresponding  vortic- 
ity patterns.  Here  again,  the  larger  transport  values  in 
the  outer  radius  are  the  result  of  the  storm  angular  momen- 
tum values  increasing  with  the  radius.  The  primary  differ- 
ence within  the  inner  volume  is  that  there  is  outward 
transport  in  the  mid  to  upper  layers  during  the  final  time 
period  in  Fig.  24b  whereas  the  transport  is  inward  in  Fig. 
19b.  This  indicates  that  the  nean  mode  contribution  to  the 
lateral  transport  is  greater  in  the  angular  momentum  budget 
than   in  the   vorticity   budget. 
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Figure  23.      Total    (solid)  ,    mean    (dashed)    and   eddy    (dash-dot) 
absolute    vorticity  lateral  transport    profiles 
for   the   0000-0600   31 T   19    February   period. 


In  summary,  the  lateral  transport  term  is  the  most 
important  term  in  the  vorticity  and  angular  momentum  budg- 
ets. The  comparisons  between  the  eddy  and  mean  mode  parti- 
tions and  the  horizontal  advection  show  that  the  cyclone's 
mean  lower  level  inflow  (convergence)  forces  the  low  tropo- 
sphere vorticity  and  angular  momentum  buildup.  At  the  upper 
levels,  the  eddy  mode  inward  transport  or  horizontal  advec- 
tion is  predominant  and  is  responsible  for  the  upper  tropos- 
pheric  increase  of  vorticity  and  angular  momentum  within  the 
vo  1  u  me  . 
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Figure  24. 


Absolute  angular  moaentum  lateral  transport  (x 
10s  cm2  sec-2)  .   Dashed  lines  indicate  negative 
values. 
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Absolute   angular   momentum   mean    mode   lateral 
transport    (x   10s   cm2    sec-2) .      Dashed   lines 
indicate   negative   values. 
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Figure  26. 


Absolute   angular   momentum   eddy    mode   lateral 
transport    (x   10s   cm2    sec-2). 
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D.   VERTICAL  REDISTR IBOTIO N 

The  vertical  redistributisn  terms  (vertical  vorticity 
advection  and  the  divergence  of  the  vertical  angular  momen- 
tum transport)  were  calculated  using  the  vertical  velocity 
fields  in  the  FGGE  data  set  (Fig.  27).  The  method  used  to 
calculate  these  values  is  briefly  described  in  the  preceding 
chapter.  Roman  (1981)  calculated  vertical  velocities  con- 
sistent with  the  lateral  transport  and  time  tendency  terms 
in  the  QLD  mass  budget  of  the  Presidents'  Day  Cyclone  (Fig. 
28)  ;  hereafter  referred  to  as  the  QLD  vertical  velocities. 
Note  that  the  QLD  vertical  velocities  are  mean  values  for 
each  layer  in  the  budget  volume  and  can  only  estimate  a  mean 
mode  vertical  redistribution  if  used  in  the  budget  calcula- 
tions. The  time  sections  of  the  FGGE  and  QLD  vertical 
velocities  reveal  significant  differences.  The  FGGE  values 
are  generally  greater  with  maximum  vertical  motions  that 
occur  earlier  and  at  lower  levels  (600mb  vice  400mb)  . 

The  FGGE  vertical  velocity  values  during  the  period  of 
most  rapid  deepening  (Fig.  29)  exceed  the  QLD  values,  espe- 
cially below  UOOmb.  The  discrepancy  between  the  FGGE  and 
mass  budget  vertical  motion  profiles  is  surprising.  If  one 
assumes  that  the  FGGE  vertical   motions  are  derived  from  the 
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Figure  27. 


FGGE   vertical   velocity  time   sections    (^bar 
sec~l) -      Dashed  lines   indicate    negative   values, 
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adiatatic  version  of  the  nonlinear  normal  mode  initializa- 
tion, the  mass  budget  omegas  might  be  expected  to  be  larger 
since  they  include  diabatic  contributions.  More  study  is 
required  to  discover  the  causes  for  these  differences  in 
vertical    motion    estimates. 

The  vertical  vorticity  advection  time  sections  (Fig.  30) 
reveal  intricate  but  similar  patterns  at  Bach  radius.  Max- 
ima occur  nearly  at  the  level  of  the  maxima  in  the  FSGE  ver- 
tical velocity  fields  (used  to  calculate  this  term)  but 
during  a  later  time  period.  Weak  downward  advection  is 
indicated  in  the  mid-layers  with  weak  upward  advection  in 
the  upper  layers.  The  mid-layer  downward  advection  must  be 
attributed  to  a  reversal  in  the  vertical  vorticity  gradient 
since  no  downward  tropospheric  vertical  /elocities  are  seen 
in  the  layer/area  averaged  FG3E  data  and  a  weak  vorticity 
decrease  with  height  is  evident  in  the  Idw  to  mid-layers  in 
the  radius  six  vorticity  time  sections  (Fig.  11).  The  maxi- 
mum value  of  the  vorticity  vertical  advection  term  (7  x 
10-* o  sec-2)  is  smaller  by  a  factor  of  two  than  the  time 
tendency  term  (13  x  10~l°  sec-2)  .  These  terms  however,  are 
considerably  smaller  than  the  maximum  value  in  the  total 
lateral        transport        term    (25   x    10~10    sec-2) .  Thus,      the 
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contribution  by  the  vertical  advection   tarm  to   the   budget 
results  is  relatively  small. 

In  the  angular  momentum  budget,  the  vsrtical  redistribu- 
tion occurs  through  the  vertical  divergence  of  the  vertical 
angular  momentum  transport  (Fig.  31) .  Significant  differ- 
ences exist  between  the  inner  and  outer  volume  patterns. 
Within  radius  six,  a  net  angular  momentum  loss  is  indicated 
in  the  two  lowest  layers,  particularly  after  the  0300  GMT  19 
February  time  period.  However,  within  the  9°  volume,  the 
pattern  reveals  a  gain  in  the  lowest  layers  and  a  loss  in 
the  450mb  layer  during  the  most  intense  period  of  surface 
development  (1200-1800  GMT  19  February).  Both  volumes  have 
mid-layer  gains  throughout  the  4  8-hour  period,  with  signifi- 
cant gains  in  the  upper  layers  in  the  inner  volume  that 
reach  a  maximum  during  the  final  period.  The  larger  volume 
pattern  shows  a  net  loss  in  the  upper  layers  throughout  the 
period  of  study  with  a  maximum  during  the  1500  GMT  19  Febru- 
ary period. 

The  angular  momentum  vertical  transport  (Fig.  32) 
through  each  level  was  calculated  to  aid  in  resolving  the 
primary  features  of  the  vertical  divergence  term.  These 
patterns  show  that   the  maximum  upward  transport   during  the 
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Figure  29. 


FGGE  (solid)  and  QLD  (dashed)  vertical  velocity 
profiles  for  1200-1800  GMT  19  February  1979. 
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initial  time  period  rises  from  850mb  to  the  700mb  level  by 
the  0900  GMT  19  February  period  in  the  6°  radius.  More 
importantly  however,  is  the  dramatic  height  increase  of  the 
zero  isopleth  during  the  early  stages  of  oyclogenesis  in  the 
inner  radius.  The  depth  of  the  upward  angular  momentum  ver- 
tical transport  reaches  a  peak  during  the  most  intense 
development  period.  Thus,  a  vertical  divergence  (loss)  of 
angular  momentum  occurs  in  the  low  troposphere  while  a  con- 
vergence (gain)  of  this  property  is  found  in  the  mid  and 
upper  troposphere. 

The  9°  pattern  is  considerably  more  detailed.  The  most 
striking  feature  is  the  downward  transport  through  all  lev- 
els during  the  1500  and  2100  GMT  19  February  periods.  This 
indicates  the  presence  of  dowaward  vertical  velocities  (not 
seen  in  the  area-average  omegas)  in  the  outer  radius  donut. 
These  downward  values  are  related  to  the  split  in  the  maxi- 
mum of  area-averaged  FGGE  vertical  velocities  (Fig.  27)  and 
the  smaller  vertical  transport  values  in  the  outer  radius 
compared  to  the  inner  radius.  The  product  of  the  downward 
velocities  and  the  larger  angular  momentum  values  at  outer 
radii  is  large  enough  to  overcome  the  6°  radius  area-aver- 
aged upward  transport  and  produce  the  downward  transport  in 
the  9°  pattern  during  this  interval. 
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In  a  general  sense,  the  relative  contribution  of  the 
vertical  redistribution  term  (vertical  divergence  of  the 
vertical  transport)  to  the  angular  momentum  budget  can  be 
seen  by  comparing  maximum  values.  The  largest  value  of  this 
term  in  the  inner  volume  (12  x  10s  cm2  see-2)  is  two  thirds 
of  the  maximum  inner  volume  time  tendency  term  (18  x  105  cm2 
sec-2).  As  in  the  vorticity  budget,  the  lateral  transport 
is  the  dominant  term  with  a  maximum  value  of  35  x  10s  cm2 
sec-2  in  the  inner  volume. 

The  analyses  in  this  sectioa  have  shown  the  upward  ver- 
tical redistribution  of  vorticity  and  angular  momentum  from 
the  low  to  the  upper  troposphere.  In  the  vorticity  budget, 
the  large  lower-tropospheric  vorticities  are  advected  upward 
although  this  term  plays  a  relatively  small  role  in  fhe 
budget.  In  the  angular  momentum  budget,  the  vertical  diver- 
gence of  the  transport  is  generally  negative  in  the  lower 
layers  and  positive  aloft  and  plays  a  more  important  role  in 
the  budget.  The  FGGE  vertical  motion  estimates  were  found 
to  differ  significantly  from  the  mass  budget  estimates. 
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Figure  31. 


Divergence  of  absolute  angular  momentum  vertical 
transport  (x  10s  cm2  sec-2) .   Dashed  lines 
indicate  negative  values  (losses). 
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Figure  32. 


Absolute  angular  momentum  vertical  transport  (x 
107  mb  cm2  sec-2) .   Dashed  indicates  negative 
values  (upward  transport). 
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E.        SOURCES/SINKS 

The  vorticity  and  angular  momentum  sinks  due  to  fric- 
tional  dissipation  were  computed  in  the  lowest  layer  using  a 
stability  independent  parameterization  scheme  (Johnson  and 
Downey,  1976)  to  ascertain  tha  significance  of  this  contri- 
bution to  the  budgets-  The  calculated  values  were  found  to 
be  negligible  compared  to  the  othar  budget  terms  and  thus 
have   no   significant  effect    on   the   budget    residuals. 

The  vorticity  tilting  term  (Fig.  33)  was  calculated  and 
its  values  are  generally  comparable  in  magnitude  to  the  time 
tendency  and  vertical  advection  values.  The  patterns  show 
weak  cyclonic  vorticity  generation  in  all  layers  during  the 
early  time  periods  and  continuing  in  the  mid-layers  through 
the      final   period.  Anticyclonic      vorticity  generation      is 

indicated  in  the  upper  and  lower  layers  during  the  1500  and 
2100  GMT  19  February  time  periods.  Strong  pattern  similari- 
ties exist  between  the  tilting  and  vertical  advection  (Fig. 
30)  terms  during  the  later  stages  of  cyclogenesis.  During 
this  time,  these  terms  exhibit  maxima  of  opposite  sign  in 
the  lower  layers  and  sign  revarsals  in  tha  mid-layers.  An 
analysis  of  the  components  of  the  tilting  term  might  reveal 
the    reasons   for   the  similaritias    and    diffarences   between   the 
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tilting  and  vertical  advection  patterns.  In  this  study  how- 
ever, the  tilting  term  was  calculated  only  to  determine  its 
significance  relative  to  the  primary  budget  terms.  In  any 
case,  the  tilting  term  patterns  are  consistent  with  the 
cyclone  study  of  Chen  and  Bosart  (1979).  Overall,  the  pat- 
terns reveal  a  net  positive  contribution  during  the  early 
stages  of  cyclogenesis  and  a  net  negative  contribution  dur- 
ing the  latter  stages. 

F.   RESIDUAL 

Residuals  were  computed  to  determine  the  integrity  of 
the  budgets.  The  vorticity  budget  residuals  (Fig.  34)  are 
near  zero  in  the  lower  layers  within  the  5°  volume  while  in 
the  9°  volume  the  values  are  relatively  small  throughout. 
The  large  negative  values  (a  vorticity  deficiency)  in  the 
upper  layers  in  the  inner  volume  during  the  6-hour  period  of 
explosive  cyclogenesis  (1200  to  1800  GMT  19  Feb)  deserve 
more  intense  scrutiny.  The  extreme  value  in  the  225mb  layer 
(approximately  -15  x  10-*°  sec~2)  is  nearly  60*  of  the  maxi- 
mum in  the  total  lateral  transport  term  (nearly  25  x  10_1° 
sec  -2)  in  the  inner  radius.  An  equally  large  positive  res- 
idual (vorticity  excess)  is  found  in  the  highest  layer  dur- 
ing   the   0300   GMT  20     February   period.         Errors   introduced   by 
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interpolation,  finite  differencing,  and  spatial  and  time 
averaging  can  contribute  to  the  residuals.  These  contribu- 
tions however,  should  be  reasonably  consistent  during  all 
time   periods   and   in   all    layers. 

Since  the  largest  negative  residuals  occur  during  the 
6-hour  period  of  most  explosive  cyclogenesis  and  deep  con- 
vection (Bosart,  1981),  it  seems  likely  that  the  intense 
mesoscale  convective  processes  are  not  adequately  resolved 
in  the  FGGE  data  set.  As  in  most  synoptic  studies,  a  major 
uncertainty  lies  in  the  vertical  velocity  fields.  The  pres- 
ence of  large  residuals  suggests  that  the  normal  mode  verti- 
cal velocities  (FGGE)  do  not  explain  the  vertical  vorticity 
redistribution   during  the    period    of  explosive   deepening. 

The  angular  momentum  budget  residual  patterns  (Fig.  35) 
are  more  complex  than  the  vorticity  residual  patterns.  The 
major  imbalances  are  again  found  during  the  period  of 
intense  development  and  convective  activity,  1200  to  1800 
GMT  19  February,  in  the  outer  volume.  The  large  negative 
values  indicate  an  angular  momentum  deficit  in  the  mid  and 
upper  layers.  Within   the  6°      volume,      a      positive   maximum 

(excess)  is  found  in  the  350mb  layer  during  the  0300  GMT  20 
February    period.      However,    the  largest   inner    radius   deficits 
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Figure  33. 


Tilting  term  in  ths  vorticity  budget  (x  10-*° 
sec-2).   Dashed  linss  indicatB  negative  value: 
Contour  interval  is  0.8  units. 


87 


are  found  in  the  lowest  layers  and  reach  a  maximum  during 
the  final  time  period.  This  indicates  a  problem  in  the 
parameterization  of  the  frictisnal  torque  for  this  data  set. 
An  interesting  point  common  to  both  budgets  is  that  the 
maximum  upper  level  deficits  occur  in  the  time  interval 
between  the  double  maxima  in  the  tendency  patterns  (Figs.  16 
and  18).  At  first  glance  this  might  indicate  a  phase  prob- 
lem in  the  budget  calculations.  Closer  examination  however, 
reveals  that  a  doubling  of  the  time  tendency  in  the  interval 
between  the  twin  maxima  is  required  to  compensate  for  the 
vorticity  and  angular  momentum  deficits  in  the  upper  levels. 
Furthermore,  this  compensation  would  require  a  substantial 
increase  in  the  magnitudes  of  the  existing  maxima  in  the 
tendency  patterns.  Thus,  it  is  unlikely  that  a  problem  in 
the  tudget  calculations  is  the  primary  cause  of  this  incon- 
sistency. Again,  the  vertical  motion  estimates  (FGGE)  used 
in  the  budget  computations  appear  suspect,  particularly 
since  the  largest  deficits  common  to  both  budgets  occur  dur- 
ing the  period  of  most  rapid  deepening. 
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Figure  34. 


Vorticity   budget   residuals    (x    10-*°    sec~2) . 
Dashed  lines   indicate   negativs    values.      Contour 
interval   is    3.0  units. 
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V.       CONCLUSIONS   AND    RECOMflENDM'IONS 

The  absolute  vorticity  and  angular  momentum  budgets  for 
the  Presidents1  Day  Cyclone,  18-19  February  1979  were  exam- 
ined employing  the  guasi-Lagrangian  diagnostic  technique  in 
pressure  coordinates.  The  FG3E  Level  Ill-b  data  set  pre- 
pared by  the  ECHWF  was  used  to  define  the  development  and 
movement    of   the    storm  as   well   as    to  calculate   the   budgets. 

The  budget  calculations  in  this  study  have  yielded  a 
number  of  interesting  results.  The  lateral  transport  is  the 
most  significant  term  in  each  budget.  Maximum  values  of  the 
lateral  transport  are  generally  twice  those  of  the  time  ten- 
dencies, which  is  usually  the  second  largsst  term.  The  mag- 
nitudes of  the  vertical  redistribution  terms  are 
approximately  half   those      of  the    time   tendency.  The    fric- 

tional  dissipation  was  small  while  the  vorticity  tilting 
term  is  equal  in  magnitude  to  the  vorticity  vertical  advec- 
tion.  The  residuals  in  each  budget  were  largest  during  the 
period  of  most  intense  convective  activity  and  explosive 
cyclogenesis  within  the  larger  volume  in  the  angular  momen- 
tum budget  and  within  the  smaller  volume  in  the  vorticity 
budget. 
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The  mass  budget  average  vertical  motion  estimates  were 
compared  with  the  area-averaged  vertical  motion  from  the 
ECMWF  analyses.  Significant  differences  were  found.  Sur- 
prisingly, the  ECMWF  vertical  motion  estimates,  which  are 
believed  to  be  from  an  adiabatic  version  of  the  nonlinear 
normal  mode  initialization,  were  larger  than  the  essentially 
kinematic   average   vertical    motions  from   the    mass   budget. 

The  mean  and  eddy  mode  partitioning  of  the  lateral 
transport  and  the  correlations  with  the  horizontal  advection 
and  divergence  components  of  the  lateral  transport  show  that 
the  mean  low-level  inflow  supplies  vorticity  and  angular 
momentum  to  the  low-troposphere.  However,  the  upper  level 
horizontal  advection  or  eddy  mode  inward  transport  is  prima- 
rily responsible  for  the  absolute  vorticity  and  storm  angu- 
lar momentum  buildup  in  the  upper  troposphere  during  the 
development   of  the  Presidents1    Day  Cyclona. 

Several  recommendations  for  further  research  emerged 
during  the  development   of  this  thesis: 

•  Vertical  velocity  fields  should  be  computed  based  on 
the  FGGE  Level  Ill-b  data  set  using  the  kinematic 
method  in  an  attempt  to  incorporate  the  diabatic  heat- 
ing contribution.  Calculation  of  the  budgets  based  on 
these  vertical  velocity  fields  could  then  be  accom- 
plished in  an  effort  to  reduce  the  residuals  and  more 
accurately   evaluate  the   budget   terms; 

•  Vertical  cross-sections  of  the  budget  terms  through  the 
storm  center  would  allow  a  more  detailed  analysis  of 
the    large-scale   dynamic   forcing   mechanisms; 


92 


The  divergence  term  of  tha  vorticity  budget  should  be 
analyzed  to  further  understand  the  relationship  between 
the  advective  and  transport  forms  of  the  lateral  trans- 
port; 

Energy      diagnostics    (kinetic     and   available      potential) 
should   also      be    completed      on   this      cyclone   for      a    more 
comprehensive   description  of   the   explosive   cyclDgenesis 
exhibited   by  this   class   of    cyclones. 
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APPENDIX    A 
COMMON    ABBREVIATIONS    AND    ACRONYMS 

cm  centimeters 

ECMWF        European  Center   for    Medium   Range    Weather 
Forecasts 

FGGE  First   GARP   Global  Experiment 

GARP  Global   Atmospheric    Research   Program 

GMT  Greenwich   Mean  Tims 

LFM  Limited   area    Fine-Mesh 

m  meters 

mb  millibars 

NVA  Negative  Vorticity    Advsction 

NMC  National   Meteorological   Center 

PVA  Positive   Vorticity    Advection 

QLD  Quasi -Lagrangian   Diagnostics 

sec  seconds 
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APPENDIX    B 

QUASI-LAGRANGIAN    BUD3ET    FRAMEWORK     AND 
GENERALIZED    BUDGET    EQUATIONS 

The  purpose  of  this  appendix  is  to  present  the  formula- 
tion of  the  storm  budget  volume  coordinate  system  and  the 
generalized  budget  equations  which  can  be  applied  to  any 
desired      atmospheric   property.  The   information     contained 

herein  is  extracted  from  Wash  (1978)  as  adopted  from  Johnson 
and    Downey    (1975A)  . 

Figure  36  depicts  the  storm  budget  volume  in  spherical 
coordinates.  The  following  notation  is  relevant  to  this 
figure: 

a  Azimuthal   coordinate 

fi  Angular  radial  coordinate 


Position  vector  from   earth's 
center 


Position  vector   from    storm   center 
to   any   point    in   budget    volume 


m,l,lc  Orthogonal  unit   vectors    along 

""  ~  meridional,    azimuthal   and 

vertical   directions   of    spherical 
coordinates    for   storm 


Wind   velocity  relative    to   earth 
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W  Horizontal  velocity  of  reference 

~°  axis 


Horizontal  velocity  at  any  point 
on  the  volume  boundary. 


Table  VI  contains  a  listing  of  the  generalized  budget 
equations.  The  following  is  an  explanation  of  the  notation 
in  these  equations: 


(0-W)        Horizontal  wind  component  relative 
~£        to  budget  volume  that  is  normal  to 
lateral  boundary 


u  Vertical  velocity  in  pressure 

coordinates 


rsin£       Radius  of  budget  volume 


Pressure 


gravity 
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Quasi 

*or^on»al 

surfaces 


Earth's 
surface 


Figure   36. 


Quasi-Langrangian  storm   budget    volume  coordinate 
system    (extracted   from  Wash,    1978). 
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TABLE    VI 

GENERALIZED    QUASI-LAGR ANGIAN    BDDGET    EQUATION    IN    PRESSURE 

COORDINATES 

(after  Wash,    1978) 


100  mb       S       2  it 
p»      /  /      /         -   fr2sinS    dad8 (-dp) 

lOOOmb    0       0         g 

where  F  is  the  volume  integral  of  the  desired  budget 
property  f. 


The  budget  equation  is 

|£   «*    LTCF)+VTlF)  +S(F) 
fit 


where  the  lateral  transport  is 

lOOmb      2ir    i  , 

LT(F)»      /  /         -(U-W)afrsinB    do(-dp) 

lOOOmb    0         ^  S  '8 

and  the  vertical  redistribution  is 

lOOmb      8       2ir    i    a  - 

VT(Fl»      /  /      /         =■  f—(.a>f)r2sinB    dadS  (-dp). 

lOOOmb    0       0         9    3P 

The  source/sink  term  is 

lOOmb       8       2tt    i    ^f       o     .     «     „    .,*  ,     ,    > 
S(F)=»      /  /      /         -  TT  r2sin6    dadS  (-dp). 

lOOOmb    0      0         9    at 
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